Low-temperature transport through a quantum dot between two superconductor leads 
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- - - We consider a quantum dot coupled to two BCS superconductors with same gap energies A. 

\l ' The transport properties are investigated by means of infinite-f/ noncrossing approximation. In 

I equilibrium density of states, Kondo effect shows up as two sharp peaks around the gap bounds. 

Application of a finite voltage bias leads these peaks to split, leaving suppressed peaks near the 
edges of energy gap of each lead. The clearest signatures of the Kondo effect in transport are three 
C~| I peaks in the nonlinear differential conductance: one around zero bias, another two at biases ±2A. 

, This result is consistent with recent experiment. We also predict that with decreasing temperature, 

the differential conductances at biases ±2A anomalously increase, while the linear conductance 
\^ ' descends. 

PACS numbers: 74.50.+r, 72.15.Qm, 73.50.Fq, 73.21.La 
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y , I. INTRODUCTION 

^ ■ In recent years, the Kondo effect for various mesoscopic systems has been extensively studied in literature. Most 
c/2 [ investigations have been concerned with a quantum dot between two normal metallic leads.^ In density of states, it has 
been found an emergence of a sharp Kondo peak near Fermi surface^ and a split of this peak when out of equilibrium .A 
In the voltage-dependent differential conductance, there appears an anomalous zero-bias maximum, ascending with 
a decrease of temperaturei^ The further studies indicated that the manifestation of the Kondo effect is strongly 
dependent on the properties of leads, by which the quantum dot is sandwiched. Since in BCS superconductors there 
are present energy gaps, it is very interesting what the Kondo effect exhibits in the case of two superconductor leads. 
Q How are transport properties modified? 
O For a quantum dot between two superconductors, recent studies have shown that the Kondo effect will cause two 

sharp peaks near gap edges in the equilibrium density of statesi^iSii This conclusion agrees with the investigation for 
dilute magnetic alloy superconductorSf^ In superconductors, all conduction electrons within energy gap are paired 
^ , into Cooper singlets and the Kondo-correlated states only can form around the energy-gap bounds. In addition, since 
00 ■ the leads coupled to a quantum dot are easily biased to nonequilibrium, it is of course raised the question, out of 
, equilibrium, how the Kondo effect affects on the density of states. Unfortunately, until now, the study about this 
' subject has not been carried out. 

i In theory, the transport properties have not been well studied yet. Ivanov predicted two surprising maxima in 
• current near biases ±2A.^ Moreover, Golub et al. only carried out the current for voltage bias in subgap regimejiS 
Very recently, the nonlinear transport through a quantum dot connected with two superconductors has been reported 
"p-, in experiment and rich phenomena were observedi^^ In the Kondo regime, except the peaks coming from Andreev 
d reflection, it has been found two peaks at biases ±2A and one peak at zero-bias in the differential conductance. The 
^ \ last one has been explained as the result of Kondo effect and the nature of the remaining two peaks stayed unclear, 
I ' however. 

■ O \ In present paper, we use the infinite-?/ noncrossing approximation (NCA) to investigate the nonequilibrium Kondo 
5 . effect for a quantum dot coupled to two superconductors. It is well known, to study the Kondo problem, NCA is a 
powerful tool and even can yield quantitative results j'^i^d^ We find that a voltage bias between two superconducting 
leads will cause the equilibrium Kondo peaks in density of states to split. The resulting peaks appear at the edges of 
energy gap of each lead and their amplitudes are suppressed. In the nonequilibrium differential conductance versus 
k> , bias, there are observed two anomalous peaks at biases ±2A and a small peak at zero bias. Our study indicates that 
5_j ■ the peaks at biases ±2A, observed in experiment, also come from the Kondo effect. 

\ This paper is organized as follows. In section II we introduce the model and NCA for the case of the leads having 
energy gaps. In section III the numerical results are discussed. Finally, the conclusions are present in section IV. 
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II. MODEL AND FORMULATION 



We consider a quantum dot coupled to two (left and right) BCS-superconducting leads. Between lead and quantum 
dot a dc voltage eVa (a=L or R) is applied. In order to clearly explore the effect of energy gaps on the Kondo 
resonance, here, we restrict oursclf to the case where the superconducting gap energy A is much larger than Kondo 
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temperature Tk- Hence, the influence of the coupHng between leads and the quantum dot on the superconductivity 
can be neglected. 

After performed a gauge transformation^^, the Hamiltonian of the system can be written as 

H= H^ + Hdot+ ^T{r)- (1) 

a=L,R a=L.,R 

Here Ha is the Hamiltonian of the superconducting lead a, 

= H'^ak^^akaCak^ + Xl(^"^"kiCa-kT + A„c+_k^c+k^), (2) 
kcr k 

with c^j,g. (c^kcr) creating (destroying) an electron with momentum k and spin a in the lead a. The dot is modelled 
as an infinite-C/ Anderson impurity having bare energy Sd- Its Hamiltonian in the slave-boson representation is given 
by 

Hdot = edYfU.- (3) 

cr 

The U oo constraint of single occupancy takes the form Q — fa fa + b'^b — 1. An auxiliary boson operator b'^ 
and fermion operator f^ create an empty site and singly occupied state, correspondingly. They relate to the ordinary 
electron creation (annihilation) operator (da) in quantum dot by 

da = b+fa, (4) 

4 = f^b. (5) 
The hopping between the lead a and the dot is described by i?^, 

(r) - J2 ^"k.e'^° c+^ab^fa + H.c. (6) 

kcr 

with (f)a = eVaT /h. 

The NCA is based on a self-consistent Feynman propagator expansion. '^•''•^^ To make the standard diagram tech- 
niques valid, for example the Wick's theorem, the slave boson perturbation expansion is initially formulated in the 
grand canonical ensemble, i.e. in the enlarged Hilbert space of pseudo-fermion and slave boson degrees of freedom. 
In a second step, the exact projection of the equations onto the physical Hilbert space, Q — 1, is performed. 

The noncrossing approximation is accurate enough to study the dynamics of the quantum dot between two normal 
leads. In this approach, only the lowest-order diagrams for the boson and fermion propagators need to be considered. 
When one or both leads are replaced by superconductors, additional diagrams should be taken into account to include 
the effect of multi-Andreev reflection.^"* In our study, we restrict ourself to the case where the Coulomb interaction 
U in quantum dot is much larger than the gap energy A and the coupling constant between the dot and leads F, 
i.e. in the infinite U limit. Due to the energetically unfavorable of the 2e-charge fluctuation, the contribution of the 
Andreev reflection can be neglectedii^ Hence, the usual noncrossing approximation also can be used here. 

The resultant equations for the self-energies of the retarded Green functions of the pseudo-fermions, G^j^{uj) — 
(ijj — £d — Syg.(ix'))~* and the slave-bosons, D^{uj) = {cu — H''(cl'))~*, constrained to the physical subspace, can be 
written as 

- ^ + ey„)(l -f{Lj-e + eVa))D-{e), (7) 



n'(^) = E / V^""^^^ " " ^^"^^^^ " " eVa)G}^{e). (8) 

aa 

Here, Vaa{'^) ~ ^aaPa\i^)- The normalized density of states of electrons in superconductors pa\uj) has the formic 

p(^)ic.)^Re^P=. (9) 
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Coupling constant Taa is defined \Vaka\Na{0) With Na{0) being the normal density of states in the Fermi 

level. 

Noted that the density of states of electrons in superconductors can be written as Eq. 0, only when the gap 
energy A is much larger than the Kondo temperature. In another case, it needs to be determined self-consistcntly, as 
shown in Ref.^ for a magnetic impurity in superconductor. 

In the nonequilibrium situation, additional Green's functions named less Green's functions of pseudo-fermions 
GJ^ {to) , and slave bosons (uj) , should be considered. Standard manipulation of the nonequilibrium Dyson equations 
leads to 

Gf^iu;) ^ G%iu;)^%{iu)G%{L,), (10) 
D<{lu) = D'-{uj)U<{uj)D''{uj), (11) 

where the self-energies are given by 

^/<x(^) = - II / —^o^aie -io + eV^)f{u - e + eV^)D<{e), (12) 

n<(c.) = -5] j '^T^,{e-oj-eV^){l-f{e-uj^eV^))G%{e). (13) 
The physical impurity Green's function is then obtained from 

GS.(c^) = | / ^[lmD^{e)G%{io + e)-D<{e)lniG}^{u^ + e)] (14) 



with the partition function 



Z^^j dc.[i?<(c.)-^G<,(.;)]. (15) 



III. NUMERICAL RESULTS 

We have numerically solved the self-consistent NCA equations. For the evaluations, it is assumed a Gaussian 
density of states of normal conduction electrons with half width D. The single level of the quantum dot is chosen 
to be £d = —0.5D. We only treat with a symmetric structure: two s-wave superconducting leads, characterized by 
A = Al = Afl = 0.25D, are coupled to a quantum dot with coupling constant F = Flo- = ^r/.a = 0.1251?. In this 
case the Kondo temperature for normal leads Tk is estimated to be 0.0171?. 

It is well known that, for a quantum dot between two metallic leads, the set of NCA equations can yield correct 
sharp peak at the chemical potential in the equilibrium density of states. This peak at zero temperature should be 
understood as the exhibition of the transition from the A^-particle ground state to the ground state with TV -|- 1 or 
N — 1 electrons!^ These iV ^ + 1 and N ^ N — 1 transitions occur only when the correlated ground state has a 
finite amplitude to have an empty site. Since by definition the ground state energies differ by the chemical potential, 
the Kondo peak in density of states should be formed at the Fermi level. 

When the leads become superconductors with A ^ Tk, all electrons within the gap are paired into Cooper singlets. 
The occurrence of the N N+1 and N N — 1 transitions should cost additional energy of order ±A. Consequently, 
in the impurity density of states, the sharp Kondo peaks are expected to be observed near the gap bounds. Indeed, 
from Fig. 1, these two peaks at ±A are clearly visible and a gap opens around Fermi surface. We also can see that 
the sharp peak at —A shifts towards the center of gap. The shape of this peak is more sharp and its amplitude is 
more large, by comparison with the peak at A. When temperature decreases, both peaks become more pronounced. 
The small weights of these two peaks refiect the fact that the probability of site being unoccupied in the ground state 
is small. Noted that with lower temperature, the shape of the peak around —A is very sharp and the numerical 
calculation is difficult to be carried out. 

In the impurity density of states, there is another low, broad peak around bare-level energy. This peak is associated 
with the transition of the empty site to excited states. For finite interaction energy ?7, an additional broad peak in 
the density of states near Sd + U should appear. At the same time, the peak around A also will move towards the 
center of gap, as that shown in Ref. SJ. 
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FIG. 1: The equilibrium density of states for a quantum dot coupled to two superconducting leads at several temperatures: 
T = IOTk, 5Tk and 4Tk. The inset shows a blowup of the region around w = —A. 
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FIG. 2: The nonequilibrium density of states for a quantum dot between two superconducting leads. Several biases are applied: 
eV/2 = 0, A/3 and A. The temperature is T = 4Tk. The inset shows a blowup of the region around u> = — A. 



We also have carried out the nonequilibrium density of states of impurity. The results are plotted in Fig. 2. It can 
be seen that when a voltage bias bVl = — cVr = eV/2 is applied, the equilibrium peaks at ±A split and the resulting 
peaks are formed around ±eV/2 ± A. In comparison with the corresponding equilibrium peaks, the amplitudes of 
these split peaks decrease and their shapes become more broadening. For the bias eV/2 = A, there only three peaks 
are observed, however. The split peaks near eV/2 — A and —eV/2 + A are superposed and a more broadening peak 
appears. Noted that, the positions eV/2 — A and —eV/2 + A also correspond to the lower edge of left-lead gap and 
upper edge of right-lead gap, respectively. 

Physically, out of equilibrium, the quantum fluctuation will produce a finite probability of an empty site. Hence, 
the Kondo peaks in density of states also can occur at the two edges of energy gap of each lead, i.e. near ±A ± eV/2. 

We also have carried out nonlinear differential conductance after computing the current by formula 

Il{V) = -JiT.J ^^L-i"" - t^LmU^) - ImG(;)(u;)/(a; - Ml)], (16) 

^«(^) = I E / '^^^.(a; - f^R)[G<^{u;) - lmG<^J^ (u) f {w - nn)], (17) 
where the lesser Green function of the impurity can be obtained from the pseudo-fermion and slave boson Green 
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FIG. 3; The voltage dependence of differential conductance for a quantum dot connected with two superconductors. The inset 
shows a blowup of the region around zero bias. 



functions via 

Gti^) = I / §^G%ie)D^e - u). (18) 

Our study is limited to the case of infinite U . The contribution of Josephson effect to current is small and can be 
neglected 1^ 

Since the NCA is a conserving approximation, the currents computed from the left and the right leads should be 
the same. In calculation, the current conservation within the NCA has been checked and the agreement between two 
currents within 0.5% is found. We evaluate the differential conductance, G{V) = dI{V)/dV , by numerical derivative. 

We plot the nonlinear differential conductance versus bias in Fig. 3 for several temperatures. A small peak at zero 
bias always can be seen. In contrast to the case of normal leads, this zero-bias maximum descends with decreasing 
the temperature. It is also observed two anomalous peaks in differential conductance at biases eV = ±2A. When 
temperature falls, the shapes of these two peaks become more sharp and their amplitudes increase. In the following, 
we will clarify that all these three peaks are produced by Kondo effect. 

To explain the absence of anomaly of zero-bias conductance, we rewrite the linear differential conductance as 

= 0) = T^^r- / du.p^i\u) (-^) l^GM. (19) 

Since p\^\lo) only has a nonzero value when out of energy gap, the contributions of the impurity spectrum and the 
Fermi- function differential to conductance come from the energy regime where > A. The lower-energy Kondo peak 
in equilibrium density of states appears within energy gap (see Fig. 1), so its contribution can be neglected. When the 
temperature falls, the higher-energy Kondo peak increases slowly and the exponential diminution of Fermi- function 
differential becomes dominant. In results, linear differential conductance is suppressed. Noted that the absence of 
anomaly of zero-bias conductance also has been carried out for a quantum dot connected with a metallic and a 
superconductor leadiii 

When the bias deviates from zero, splitting of equilibrium Kondo peaks in density of states and suppression of 
the resultant peaks in amplitudes should cause a decrease of the differential conductance. However, for \eV\ « 2A, 
from Fig. 2 we can see, in the density of states, there appears a peak produced by superposition of the higher-energy 
and lower-energy Kondo peaks. The position of this peak coincides with the edges of left and right leads, so perfect 
transparency occurs. Hence, in differential conductance, two peaks around biases ±2A emerge. Furthermore, with 
decreasing temperature, the enhancement of the amplitude of this superposition peak leads to a rapid increase of 
nonlinear conductance near biases ±2A. 

In the experiment, besides the peaks induced by Andreev reflection, such three peaks in voltage dependence of 
nonlinear conductance also have been observed^ However, there, the origin of the peaks at ±2A was not clear. 
Our calculation indicated these two peaks are produced by Kondo effect. Unfortunately, the quantitative comparison 
between our results and experiment can not be made, since our study for A 3> Tk corresponds to the experimental 
situations, where have two energy levels of quantum dot down to Fermi surface. 
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IV. CONCLUSIONS 



We have studied the nonequihbrium Kondo effect for a quantum dot coupled to two superconductors with energy 
gap A ^ Tk- Employing the infinite-C/ noncrossing approximation, we have carried out the nonequihbrium impurity 
density of states. When a voltage bias is applied to the both side of leads, the peaks, produced by Kondo effect 
at equilibrium, split and give rise to suppressed peaks around the bounds of energy gap of each lead. We also 
have evaluated the differential conductance. The Kondo effect in the voltage dependence of differential conductance 
shows up as two anomalous peaks around biases ±2A and a small peak near zero bias. This result accords with the 
experimental observation. With descending the temperature, peaks at biases ±2A ascend and the zero-bias peak falls. 
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